Multichannel e ± scattering on excited Ps states 
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Scattering and reaction cross sections of e — Ps system are calculated for total angular momentum 
L = 0, 1 and 2 and energies between the Ps(n = 2) — Ps(n — 3) threshold. We solved a set of 
Faddeev-Merkuriev and Lippmann-Schwinger integral equations by applying the Coulomb-Sturmian 
separable expansion technique. We found that the excited positronium states play dominating roles 
in scattering processes. 
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The e — Ps system plays a very important role in 
studying the antimatter. While, on the experimental 
side, new experiments, that involve the positronium in 
one way or another, are being carried out or planned 
jjj, on the theoretical side, the existing calculations are 
restricted for low energy elastic scattering below the 
Ps(n — 2) threshold (see Ref. || and references therein). 

Recently, based on a three-potential picture, we have 
developed a new method for treating three-body Coulom- 
bic systems ||. The three-potential formalism results in 
a set of Faddeev-Merkuriev and Lippmann-Schwinger in- 
tegral equations. These integral equations were solved by 
the Coulomb-Sturmian separable expansion technique. 

In this paper we present scattering and reaction calcu- 
lations for energies between the Ps(n = 2) — Ps(n = 
3) threshold and for total angular momentum L = 
0,1, and 2. First we outline the method of Ref. || to 
the e — Ps system and then present the results. 

In the e ± — Ps system two particles are always identi- 
cal. Let us denote them by 1 and 2, and the non-identical 
one by 3. The Hamiltonian is given by 



H = PT 



"3 i 



(1) 



where H is the three-body kinetic energy operator and 
!)„ denotes the Coulomb interaction in the subsystem a. 
We use the usual configuration-space Jacobi coordinates 
x a and y a ; x a is the coordinate between the pair (/?, 7) 
and y a is the coordinate between the particle a and the 
center of mass of the pair (/?, 7). Thus the potential u„, 
the interaction of the pair ((3, 7), appears as v„ (x a )- We 
also use the notation X = {x a ,y a } E R 6 . 

The Hamiltonian ([!]) is defined in the three-body 
Hilbert space. So, the two-body potential operators are 
formally embedded in the three-body Hilbert space, 



! (x)ly, 



(2) 



where l y is a unit operator in the two-body Hilbert 
space associated with the y coordinate. The role of a 
Coulomb potential in a three-body Coulombic system 
is twofold. In one hand, it acts like a long-range po- 
tential since it modifies the asymptotic motion. On the 



other hand, however, it acts like a short-range potential, 
since it correlates strongly the particles and may support 
bound states. Merkuriev introduced a separation of the 
three-body configuration space into different asymptotic 
regions Q . The two-body asymptotic region Q is defined 
as a part of the three-body configuration space where the 
conditions 



\x\ <x (l + M/ 2 /o) lA ', 



(3) 



with xo,yo > and v > 2, are satisfied. Merkuriev 
proposed to split the Coulomb interaction in the three- 
body configuration space into short-range and long-range 
terms 



v c = 



,(0 



(4) 



where the superscripts s and / indicate the short- and 
long-range attributes, respectively. The splitting is car- 
ried out with the help of a splitting function £, 



v (s \x,y) = v C (x)((x,y), 
v^(x,y) = v c (x)[l-ax,y)}. 

The function £ is defined such that 



x- 



1, x en 

otherwise. 



(5) 
(6) 



(7) 



In practice usually the functional form 

C{x, V) = 2/ {1 + exp [(x/x°r/(l + y/y )] } , (8) 



is used. 

In the Hamiltonian ([!]) the potential v§ , acting be- 
tween the identical particles, is a repulsive Coulomb po- 
tential which does not support bound states. Conse- 
quently, the entire v§ can be considered as long-range 
potential. Then, the long-range Hamiltonian is defined 
as 



,(0 , „,(0 



(9) 



1 



and the three-body Hamiltonian takes the form 
H = H®+vP+vi'\ 



(10) 



So, the Hamiltonian of the — Ps system appears for- 
mally as a three-body Hamiltonian with two short-range 
potentials. The Faddeev procedure is applicable, and, 
in this case, we get a set of two-component Faddeev- 
Merkuriev integral equations 



i^} = i$f ) ) + GN s V 2 } 
IV> 2 ) = <?<'>«<"> Ilk), 



(11) 

(12) 



where \<pa^} and Gq' are eigenstate and resolvent opera- 
tor, respectively, of the channel Coulomb Hamiltonian 



(13) 



Since the particles 1 and 2 are identical the Faddeev 
components \ip\) and {tfa), in their own natural Jacobi 
coordinates, have the same functional forms 



Therefore we can determine \ip) from the first equation 
only 



^) = \^ ) ) + G^v[ s >pV\^), 



(').,(-). 



(15) 



where V is the operator for the permutation of indexes 
1 and 2 and p = ±1 are eigenvalues of V. We note that 
although this integral equation has only one component 
yet gives full account on the asymptotic and symmetry 
properties of the system. 

We solve this integral equation by using the Coulomb- 
Sturmian separable expansion approach. The Coulomb- 
Sturmian (CS) functions are defined by 



(r\nl) 



-1 1/2 



(n + 2l+l)) 



(2br) l+1 exp(-br)L% +1 (2br), 

(16) 



with n and I being the radial and orbital angular mo- 
mentum quantum numbers, respectively, and b is the 
size parameter of the basis. The CS functions 
form a biorthonormal discrete basis in the radial two- 
body Hilbert space; the biorthogonal partner defined by 
(r\nl) = (r\nl)/r. 

Since the three-body Hilbert space is a direct product 
of two-body Hilbert spaces an appropriate basis can be 
defined as the angular momentum coupled direct product 
of the two-body bases 

\nvl\) a = \nl) a ® \vX) a , (n, v = 0, 1, 2, . . .), (17) 

where \nl) a and |i / A) Q , are associated with the coordinates 
x a and y a , respectively. With this basis the complete- 
ness relation takes the form (with angular momentum 



summation implicitly included) 

\nvlX) a a (ni>lX\ = lim 1 



N 

lim > 

n : u— 



N 



(18) 



where (xy\nvlX) = (xy\nvlX) / (xy) . 

We make the following approximation on the integral 
equation (|l5|) 



\i,)=\^)+GT^vr P vim 



(s) 

i.e. the operator v\ pV is approximated in the three- 
body Hilbert space by a separable form, viz. 

v[ s) pV = lim l?v[ s) pVl? 

N^oo 



(19) 



l?v[ s) pVl? 



N 



MA)i«i ,) i(fiVI / A'| ) (20) 

= i{nvlX\iff ) pV\n'v'l'X')i. Utilizing the 
properties of the exchange operator V these matrix 



(14) where 



elements can be written in the form 



,(») _ 



p x 



i{nvlX\v^ \n'v'l'X')2, and can be evaluated numerically 
by using the transformation of the Jacobi coordinates. 

Now, with this approximation, the solution of the in- 
homogeneous Faddeev-Merkuriev equation turns into a 
solution of a matrix equation for the component vector 
ip = 1 (nvlX\tp) 



C0> 



where 



and 



dp = iinvlXlG^ln'v'l'X')!. 
The formal solution of Eq. ( pil| ) is given by 
[(G i °)- 1 -^% = (Gf)- 1 $«. 



(0 



(nvlX\$f) 



3(0 I 



(21) 



(22) 



(23) 



(24) 



Unfortunately neither Gjp nor are known. They 



(0 



are related to the Hamiltonian , which itself is a com- 
plicated three-body Coulomb Hamiltonian. In the three- 
potential formalism H G^ is linked to simpler quantities 
via solution of a Lippmann-Schwinger equation, 

(Gjty^&r 1 -!/}, (25) 

where 

Si^w^x' =iMA|G 1 |nWA')i (26) 



2 



and 



U 1 

— ni/lX.n u I A 



i{nvl\\U r \rii/l'\')-i 



(27) 



The operator G\ is the resolvent operator of the Hamil- 
tonian 



Hi =H°+v?. 
The polarization potential U 1 is defined by 



U 1 = v. 



(i) 



u 3 > 



(28) 



(29) 



and its matrix elements can again be evaluated numeri- 
cally. 

Similarly, also $j can be linked to simpler quantities 

[(Si)- 1 -E 1 ]ti ,) = (Si)- 1 Ii, (30) 

where ^ = i (nz/i A |<I>i), and 4>i is eigenstate of Hi. 

The three-particle free Hamiltonian can be written as 
a sum of two-particle free Hamiltonians 



hi 



(31) 



Consequently the Hamiltonian Hi of Eq. (pq) appears as 
a sum of two Hamiltonians acting on different coordinates 



H x = h Xl 



"VI ) 



(32) 



with h Xl = h xi +Vi(xi) and h Vl = h yi , which, of course, 
commute. Therefore its eigenstate, in CS representation, 
appears as 



y(nvl\\$i) = i(nl\(j>i) x i{v\\xi), 



(33) 



where \<pi) and |xi) are bound and scattering eigenstates 
of h Xl and h yi , respectively. 

The matrix elements of Gi can be determined by mak- 
ing use of the convolution theorem 



Gi(z) 



1 
2vri 



dz'g xi (z~z')gjz'), (34) 



where g xi and g yi are resolvent operators of h xi and h yi , 
respectively. The corresponding CS matrix elements of 
the two-body Green's operators for all complex energies 
and of the two-body solutions in Eq. (|3^) are known an- 
alytically. The contour C should encircle, in positive di- 
rection, the spectrum of h Vl without penetrating into the 
spectrum of h Xl . Further details on the contour and on 
those CS matrix elements are given in Ref. || and refer- 
ences therein. 

The S'-matrix of the — Ps scattering process, in the 
three-potential picture |3J, can be decomposed as 



S 



(2) 
fi 



-2mS(E f - Ei) 

x({$t ) |C/ 1 |^ )(+) ) + (^ ( " ) l^ ) l4! ) }), (35) 



where i and / refer to the initial and the final states, 
respectively. Having the solutions ip and $^ the matrix 
elements can easily be evaluated. 

In the numerical calculations we used atomic units. 
For the parameters of the splitting function (||) we took 
xq = 10, yo = 20 and v = 2.1, respectively, and for the 
size parameter of the CS basis we used b — 0.2. In the 
expansion of the potentials we went up to 11 angular mo- 
mentum channels and, in each angular momentum chan- 
nels, up to N = 27 CS functions. This way we achieved 
convergence up to 2 — 3% and also the iT-matrix were 
symmetric with a similar accuracy. Some of the results 
were cross-checked by the results of configuration-space 
differential equation calculation 0, and we found again 
very good agreements. 

The results for total angular momentum L = 0, L = 1 
and L = 2 are given in Tables @, @ and |h| respec- 
tively. We can see that the excited positronium states 
play dominating roles in scattering processes, especially 
when the total energies approach the positronium exci- 
tation threshold (from above). This behavior is consis- 
tent with the rather large size of the excited positron- 
ium targets, where the long-range polarization potential 
play dominant roles. This behavior is similar to the phe- 
nomena found in p — Ps multichannel scattering process, 
where this mechanism dominates the antihydrogen for- 
mation cross section (will be published in a separate pa- 
per). 
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TABLE I. L = partial cross sections (in na^) in the 
Ps(n = 2) — Ps(n = 3) gap. Channel numbers 1,2 and 3 

0) +Ps(2s) 

Ps(2p) respectively. 



denote the channels e ± (A = 0) + Ps(ls), e (A 
and e ± (A = '. 



fcl 


Ch.# 


1 


2 


3 


L = p = l 




1 


4.384 


0.043 


0.032 


0.51 


2 


1.125 


67.79 


0.245 




3 


0.851 


0.247 


301.6 




1 


4.055 


0.043 


0.033 


0.52 


2 


0.576 


60.77 


45.45 




3 


0.445 


45.51 


62.48 




1 


3.756 


0.043 


0.037 


0.53 


2 


0.402 


43.57 


19.28 




3 


0.349 


19.32 


63.77 




1 


3.484 


0.044 


0.041 


0.54 


2 


0.324 


28.39 


2.12 




3 


0.286 


2.13 


71.74 


L = p=-l 




1 


12.52 


0.0001 


0.0000 


0.51 


2 


0.0023 


121.69 


0.941 



3 







3 


0.0011 


0.936 


263.84 






1 


11.79 


0.0001 


0.0001 


0.52 




2 


0.0021 


102.42 


35.32 






3 


0.0009 


35.28 


42.75 






1 


11.10 


0.0002 


0.0002 


U.Oo 




Z 


U.UUzo 


OO.OU 








3 


0.0019 


31.28 


37.06 






1 


10.45 


0.0003 


0.0004 


0.54 




2 


0.0020 


8.93 


19.78 






3 


0.0025 


19.75 


44.23 


TABLE II. 


L = 


1 partial 


cross sections (in ttOq 


) in the 


Ps(n 


= 2)- 


Ps(n 


= 3) gap. 


Channel numbers 1,2,3 and 4 


denote the channels 


e ± (A = l) 


+ Ps(ls) 


,e±(A = l) + 


Ps(2s), 


e ± (A 


= 0) + Ps(2p) and e ± (A 


— 2) + Ps(2p), respectively. 


fei 


Ch.# 


1 


2 


3 


4 








L = 1 


p=\ 








1 




20.22 


0.090 


0.466 


0.246 


0.51 


2 




2.378 


296.76 


30.16 


29.21 




3 




12.21 


30.17 


90.13 


139.14 




4 




6.49 


29.40 


138.90 


591.54 




1 




19.24 


0.086 


0.617 


0.317 


0.52 


2 




1.16 


17.76 


56.03 


26.56 




3 




8.33 


55.91 


127.26 


66.25 




4 




4.28 


26.37 


66.05 


305.29 




1 




18.32 


0.095 


0.804 


0.398 


0.53 


2 




0.874 


65.19 


50.25 


38.59 




3 




7.47 


50.16 


69.26 


28.60 




4 




o.oy 


OO.OO 


ZO.OZ 


1 I 4.0O 




1 




16.91 


0.207 


0.916 


0.599 


0.54 


2 




1.48 


81.55 


50.21 


19.94 




3 




6.50 


50.12 


66.48 


10.43 




4 




4.20 


19.75 


10.36 


169.84 








L = 1 


p=-l 








1 




12.80 


0.458 


0.184 


0.175 


0.51 


2 




12.16 


290.22 


44.19 


9.62 




3 




4.89 


44.28 


737.30 


135.86 




4 




4.66 


9.63 


135.60 


813.33 




1 




11.93 


0.484 


0.201 


0.194 


0.52 


2 




6.60 


21.53 


9.78 


66.71 




3 




2.71 


9.77 


246.08 


46.32 




4 




2.64 


66.40 


46.17 


308.15 




1 




11.14 


0.514 


0.229 


0.214 


0.53 


2 




4.82 


81.03 


9.85 


79.18 




3 




2.13 


9.87 


52.74 


14.22 




4 




2.01 


78.89 


14.16 


94.31 




1 




10.41 


0.542 


0.262 


0.232 


0.54 


2 




3.96 


145.62 


4.51 


50.12 




3 




1.91 


4.52 


3.24 


4.66 




4 




1.68 


49.97 


4.65 


25.25 



TABLE III. L — 2 partial cross sections (in 7rao) in the 
Ps(n = 2) — Ps(n — 3) gap. Channel numbers 1,2,3 and 4 
denote the channels e ± (A = 2) +Ps(ls), e ± (A = 2) +Ps(2s), 
e ± (A = 1) + Ps(2p) and e ± (A = 3) + Ps(2p), respectively. 



fcl 


Ch.# 


1 


2 


3 


4 






L = 1 


p= 1 








1 


7.78 


0.010 


0.005 


0.004 


0.51 


2 


0.236 


917.53 


174.78 


64.21 




3 


0.113 


175.24 


873.28 


203.44 




4 


0.097 


63.72 


203.03 


888.39 




1 


7.58 


0.025 


0.004 


0.012 


0.52 


2 


0.333 


227.40 


28.16 


78.54 




3 


0.046 


28.16 


773.01 


69.57 




4 


0.153 


76.55 


69.05 


428.25 




1 


7.38 


0.048 


0.003 


0.025 


0.53 


2 


0.438 


68.33 


10.13 


92.31 




3 


0.026 


10.02 


502.06 


30.49 




4 


0.217 


89.59 


30.02 


231.33 




1 


7.18 


0.081 


0.003 


0.046 


0.54 


2 


0.558 


40.57 


6.30 


90.67 




3 


0.028 


6.29 


307.42 


15.42 




4 


0.311 


88.17 


15.02 


127.53 






L = 1 


p= -1 








1 


6.61 


0.025 


0.037 


0.026 


0.51 


2 


0.644 


941.82 


194.33 


36.50 




3 


0.950 


196.67 


79.94 


58.60 




4 


0.656 


35.57 


57.09 


1080.61 




1 


6.40 


0.040 


0.092 


0.054 


0.52 


2 


0.546 


427.67 


3.46 


16.52 




3 


1.22 


3.29 


293.60 


0.754 




4 


0.728 


15.06 


0.860 


510.53 




1 


6.18 


0.050 


0.135 


0.078 


0.53 


2 


0.465 


159.41 


37.79 


20.14 




3 


1.22 


37.34 


325.85 


7.27 




4 


0.711 


18.57 


7.22 


314.93 




1 


5.95 


0.055 


0.221 


0.112 


0.54 


2 


0.386 


72.91 


49.02 


22.41 




3 


1.55 


48.23 


256.51 


11.45 




4 


0.757 


21.33 


11.44 


210.85 
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